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Bismuth-metal graphite (MG) has a unique layered structure where Bi nanoparticles are encap- 
sulated between adjacent sheets of nanographites. The superconductivity below Tc (= 2.48 K) is 
due to Bi nanoparticles. The Curie-like susceptibility below 30 K is due to conduction electrons 
localized near zigzag edges of nanographites. A magnetic-field induced transition from metallic to 
semiconductor-like phase is observed in the in-plane resistivity pa around He (~ 25 kOe) for both 
_H"_Lc and H\\c (c: c axis). A negative magnetoresistance in pa for _ff_Lc (0< H <3.5 kOe) and 
a logarithmic divergence in pa with decreasing temperature for H\\c {H > 40 kOe) suggest the 
occurrence of two-dimensional weak localization efi'ect. 

PACS numbers: 71.24.-|-q, 74.80.Dm, 72.15.Rn, 71.30.+h 



I. INTRODUCTION 

A weak localization theory predicts a logarithmic di- 
vergence of the resistivity in the two-dimensional (2D) 
electron systems as the temperature (T) is loweredii In 
high-mobility Si metal oxide-semiconductor field-effect 
transistor (MOSFET), the in-plane resistivity for a sys- 
tem with an electron density n larger than a critical elec- 
tron density ric decreases with decreasing T , indicating 
a metallic behavior fM'^ This metallic state is completely 
destroyed by the application of an external magnetic field 
{H) applied in the basal plane when H is higher than a 
threshold field He- Such coplanar fields only polarize the 
spins of the electrons, indicating that the spin state is 
significant to the high conductivity of the metallic state. 
The scaling relation of the in-plane resistivity collapses 
into two distinct branches above and below He- Such be- 
haviors are very similar to those observed in amorphous 
ultrathin metal films of InOa;;^ MoGejS^ and Bij^ which 
undergo magnetic field-induced transitions from super- 
conducting phase to insulating phase. 

Bi-metal graphite (MG) constitutes a novel class of 
materials having unique layered structures. This system 
can be prepared from the reduction by Li-diphenylide 
from an acceptor-type BiCls graphite intercalation com- 
pound (GIG) as a precursor material. Ideally, the stag- 
ing structure of Bi-MG would be the same as that of 
BiCla GIG Aifliii In Bi-MG, Bi atoms would form inter- 
calate layers sandwiched between adjacent graphite lay- 
ers. For the stage-n (= 1, 2,- • ■) structure, the pack- 
ages of n graphite layers and one Bi layer would be pe- 
riodically stacked along the c axis perpendicular to the 
graphite basal plane. In BiGls GIG, the ratio of c-axis 
resistivity pc to the in-plane resistivity pa (A — pd Pa) 
at room temperature is (7.4 - 4.9) x 10'^,^"'^ suggesting 
that virtually all 7r-electrons are directionally localized; 
i.e., they can move freely along the graphite basal planes. 



but are unable to diffuse across the stack of layersii? Such 
a large ratio A is due to the insulating BiGls layer sand- 
wiched between the adjacent graphite layers. There is 
no overlapping of the wave functions over nearest neigh- 
bor graphite layers. The situation may not drastically 
change in Bi-MG in spite of the fact that BiGla layer is 
replaced by metallic Bi layer. In fact the c-axis resistivity 
of Bi-MG (= 0.1 ricm at 298 K) is almost the same as 
that of BiGla GIG at the same T, suggesting that Bi-MG 
behaves like a quasi 2D conductor. 

In this paper we have undertaken an extensive study 
on the transport and magnetic properties of Bi-MG. We 
show that this compound undergoes a superconducting 
transition at = 2.48 K. A magnetic-field induced tran- 
sition from metallic to semiconductor-like phase is ob- 
served in Pa around H ^ 25 kOe for HJ-c and H\\c (c: c 
axis). These results of Bi-MG are compared with those 
of Si MOSFET, and amorphous ultrathin metal films of 
InO:rr MoGe^ and Bi.^ 

Structural studies of Bi-MG reveal that Bi layers 
are formed of Bi nanoparticles which are encapsulated 
between adjacent sheets of nanographites. The size 
of nanographites in Bi-MG is much smaller than the 
in-plane coherence size of the graphite layers of pris- 
tine graphite. The superconductivity is mainly due 
to Bi nanoparticles, while the magnetism is mainly 
due to nanographites. Nanographites are nanometer- 
sized graphite fragments which represent a new class 
of mesoscopic system intermediate between aromatic 
molecules and extended graphene sheets. Fujita and co- 
workersi^ii^ii^ have theoretically suggested that the elec- 
tronic structures of finite-size graphene sheets depend 
crucially on the shape of their edges. The graphene 
edge of an arbitrary shape consists of two- types of edges, 
zigzag type and armchair type. The former has a trans- 
polyacetylene type structure and the latter has a cis- 
polyacetylene one. Finite graphite systems having zigzag 
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edges exhibit a special edge state. The corresponding en- 
ergy bands are almost flat at the Fermi energy, thereby 
giving a sharp peak in the density of states (DOS) at 
the Fermi energy. This is in contrast to the case of 2D 
graphene sheet with infinite size, where the DOS is zero 
at the Fermi energy. According to Wakabayashi ct al.,^^ 
the magnetism of nanographites is characterized by Pauli 
paramagnetism and orbital diamagnetism from conduc- 
tion electrons. They have predicted that the Pauli para- 
magnetic susceptibility for the nanographites with zigzag 
edges shows a Curie-like behavior at low temperatures, 
which is in contrast to a T-indcpcndcnt Pauli paramag- 
netism in normal metals. In this paper we show that the 
susceptibility of Bi-MG has a Curie-like behavior at low 
T. This behavior is discussed in the light of the prediction 
by Wakabayashi et ali^ 



II. EXPERIMENTAL PROCEDURE 

BiCls GIC samples as a precursor material, were pre- 
pared by heating a mixture of highly oriented pyrolytic 
graphite (HOPG) [grade ZYA from Advanced Ceramics, 
Ohio] and an excess amount of BiCla at 200 °C in a am- 
poule filled with chlorine gas at a pressure of 375 Torri2ii£ 
The reaction was continued for three days. It was con- 
firmed from (DDL) x ray diffraction (Rigaku RINT 2000 
x-ray powder diffractometer) that the BiCla GIC sam- 
ple consists of stage-2 as majority phase and stage-3 and 
stage-4 as minority phase. The c axis repeat distance is 
13.17 ± 0.05 1 for stage-2, 15.85 ± 0.25 1 for stage-3, 
and 20.22 ± 0.25 A for stage-4, respectively. No Bragg 
reflection from the pristine graphite is observed. 

The synthesis of Bi-MG was made by the reduction by 
Li-diphenylide from BiCla GIC. BiCla GIC samples were 
kept for three days in a solution of lithium diphenylide 
in tetrahydrofuran (THF) at room temperature. Then 
the samples were filtered, rinsed by THF, and dried in 
air. Finally the samples were annealed at 260 °C in a hy- 
drogen gas atmosphere for one day. The structure of Bi- 
MG thus obtained was examined by (OOL) x-ray diffrac- 
tion, and bright field images and selected-area electron 
diffraction (SAED) (Hitachi H-800 transmission electron 
microscope) operated at 200 kV. The same methods for 
the structural analysis were used for Pd-MGiiSiiLiS The 
(OOL) x-ray diffraction pattern of Bi-MG is much more 
complicated than that of BiCla GIC, which makes it dif- 
ficult to calculate the average particle thickness from the 
identity period in Bi-MG. Note that graphite reflections 
appear in Bi-MG, suggesting that a part of Bi atoms 
leaves from the graphite galleries occupied by BiCla inter- 
calate layers in BiCls GIC during the reduction process. 
Such Bi atoms tend to form multilayered Bi nanoparti- 
cles in the graphite galleries in Bi-MG. The number of 
Bi layers in possible multilayered structures could not 
be exactly determined at present. SAED pattern of Bi- 
MG consists of polycrystalline diffraction rings, suggest- 
ing that there are at least four Bi layers in thickness. Re- 



TABLE I: Experimentally observed spacings dexp for Bi-MG 
and dpDF given by powder diffraction file (PDF) 5-0519 for 
rhombohedral Bi [spacegroup: R-3m (ao = 4.546 A, co = 
11.860 A for the notation of the hexagonal close-packed struc- 
ture)]. All reflections can be indexed as (hkl) reflections of 
Bi. Graphite reflections are not included. 
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flections from Bi and graphite were observed. As listed in 
Tabled all the Bi reflections were labeled and attributed 
to the formation of rhombohedral Bi, according to the 
standard ICDD PDF (Card No. 05-0519). This result 
indicates that Bi nanoparticles are crystallized as rhom- 
bohedral Bi phase in Bi-MG. The observed spacings of 
Bi-MG are 1-2 % shorter than those of bulk Bi metal. 
The size of Bi nanoparticles distribute widely around the 
average size 110 A. More than 50 % of Bi nanoparti- 
cles has sizes ranging between 10 and 50 A. The largest 
particle size is 750 A. 

The measurements of DC magnetization and resistiv- 
ity were made using a SQUID magnetometer (Quantum 
Design MPMS XL-5) with an ultra low field capability 
and an external device control mode. In the present work 
we used a Bi-MG sample based on HOPG which is par- 
tially exfoliated. The stoichiometry of C and Bi was not 
determined. The in-plane resistivity pa and the c axis re- 
sistivity pc were measured by a conventional four-probe 
method. The sample had a rectangular form with a base 
6.0 X 1.6 mm^ and a height 0.47 mm along the c axis. 
For the measurement of pa, four thin gold wires (25-/im 
diameter) that were used as as the current and voltage 
probes were attached to one side of the c surfaces by sil- 
ver paste (4922N, du Pont). For the measurement of pc, 
two thin gold wires as the current and voltage probes 
were attached to each c-surface of the sample. The cur- 
rent (/ = 10 mA for pa and 3 mA for pc) was supplied 
through the current probes by a Keithley type 224, pro- 
grammable DC current source. The voltage V generated 
across the voltage probes was measured by a Keithley 
182 nanovoltmeter. The linearity of I-V characteristics 
was confirmed for the measurements of pa and pc. 



III. RESULT AND DISCUSSION 

A. Meissner effect due to Bi nanoparticles 

We measured the DC magnetization of Bi-MG. After 
the sample was cooled from 298 to 1.9 K at 77 = 0, the 
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FIG. 1: T dependence of ZFC and FC susceptibilities of Bi- 
MG based on HOPG. H = 100 Oe. HLc (c: c axis). The 
inset shows the magnetic phase diagram of _ff vs T for HLc, 
where the peak temperatures of ZFC and FC susceptibilities 
are plotted as a function of H. The solid line denotes the 
least-squares fit of the ZFC data (•) to Eq.Q. 



measurement of zero field cooled (ZFC) magnetization 
MzFC "was made with increasing T from T9 to 15 K in 
the presence of H. The sample was kept for 20 min- 
utes at 50 K. Then the sample was cooled from 50 to 15 
K. The measurement of field-cooled (FC) magnetization 
Mpc was made with decreasing T from 15 to 1.9 K in 
the presence of the same H . For convenience, hereafter, 
the direction of H in the c plane is denoted as i?_Lc (c: 
c axis) . Figures ^ and [5] show typical examples of the 
T dependence of the susceptibility Xzfc {=Mzfc/H) 
and xfc {=Mfc/H) at H = 100 Oe for H±c and H\\c, 
respectively. Both xfc and xzfc show a sharp peak 
around 2.5 K at = 20 Oe, which results from the 
competition between the diamagnetic susceptibility due 
to the Meissner effect and the Curie-like susceptibility 
(which will be described later). The peak shifts to the 
low T side with increasing H. No peak is observed above 
1.9 K for H > 250 Oe. For convenience, the peak tem- 
perature is defined as a critical temperature Tc{H). The 
peak temperature Tc{H) for xfc is almost the same as 
that for Xzfc- In the insets of Figs, ^andl^l we show 
the H-T magnetic phase diagram for Bi-MG for iJ_Lc 
and H\\c, respectively. We find that the data of if vs T 
(ZFC susceptibility) for ii_Lc and H\\c are well described 

by 



H = Ho{ 



Tc-T 



(1) 



Bi-MG 




FIG. 2: T dependence of ZFC and FC susceptibilities of Bi- 
MG. H = 100 Oe. H\\c (c: c axis). The inset shows the 
magnetic phase diagram oi H ys T for H^c. The solid line 
denotes the least-squares fitting curve for the ZFC data (•) 
to Eq.im. 



where is the critical field at T = K {Hq = Hc2 for 
H^-C and Tio = Hcs, for i?||c), (3 is an exponent, and Tc 
is a critical temperature at H = Q. The field Hc^ is a 
field to nucleate a small superconducting region near the 
sample surface. The least squares fits of the data of H 
vs T (ZFC susceptibility) yield the values of iJc2 = 489.4 
± 0.5 Oe, P = 0.48 ± 0.02, = 2.48 ± 0.06 K for i7±c, 
and i7c3 = 838 ± 1 Oe, /3 = 0.69 ± 0.02, = 2.48 ± 0.02 
K for H\\c. Here we note that the superconductivity is 
observed in ultrathin films of amorphous Bi grown in top 
of layer of amorphous Gc (thickness, 6-10 A). Haviland 
et a\t^ have reported that the superconductivity occurs 
for the thickness of Bi films, d = 6.73 - 74.27 °A. The 
critical temperature Tc decreases with decreasing d: Tc 
= 5.6 K for d = 74.27 A and Tc « 0.8 K for d « 6.73 
A. Markovic et ali^ have shown that Tc decreases with 
decreasing d for d > dc {dc = 12.2 A): Tc = 0.5 K for d = 
15 A. Weitzel and Micklitz^" have reported surface su- 
perconductivity in granular films built from well-defined 
rhombohedral Bi clusters (mean size = 38 A) embedded 
in difference matrices (Kr, Xe, Ge) or with H2 or O2 gas 
adsorbed on the cluster surface. The critical tempera- 
ture Tc, which is dependent on matrices and gases used, 
is between 2 and 6 K. This value of Tc is on the same 
order as that of Bi-MG. The exponent f3 for H±c is very 
close to that (/3 = 0.5) predicted for homogeneous system 
of isolated superconducting grainsiSi The ratio Hcz/Hc2 
is calculated as 1.71, which is very close to the predicted 
value 1.695.^^ Similar behavior is observed in the critical 
fields for HLc and i7||c in stage- 1 K GIG (KCg)^ The 
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FIG. 3: Hysteresis loop of DC magnetization for Bi-MG. 
HA-C. T — 1.9 K. The measurement was made with increasing 
H from to 500 Oe (denoted by •), with decreasing H from 
500 Oe to -500 Oe (O), and with increasing H from -500 Oe 
to 500 Oe (A). 



coherence length ^ is estimated as ^ = 820 A from the 
value of Hc2 using the relation Hc2 = ^q/{2tt^'^), where 
3>o (= 2.0678 X 10""^ Gauss cm^) is a fluxoid. The co- 
herence length ^ is much larger than the size of islands 
(110 A in average). 

Figure O shows the hysteresis loop of the magnetiza- 
tion Ma{H) for H±c at T = 1.9 K. The sample was 
cooled from 298 K to 1.9 K at = 0. The magnetiza- 
tion Ma{H) was measured with increasing H from to 
500 Oe, with decreasing H from 500 Oe to -500 Oe, and 
with increasing H from -500 Oe to 500 Oe. The magne- 
tization consists of superconductivity contribution and 
diamagnetic background. In Fig. QJa) we show the hys- 
teresis loop of AMa{H) that is defined by Ma{H) at 1.9 
K minus a diamagnetic background given by XdH with 
Xd = -5.143 X 10"'' emu/gOe. A hysteresis loop char- 
acteristic to a type-II superconductor is observed with a 
lower critical field Hd (= 80 Oe). In contrast, the mag- 
netization hysteresis \oopMa{H) at T = 3.3 K is very 
different from that at 1.9 K. It seems that there is nei- 
ther local minimum nor local maximum. In Fig. 2fb) 
we show the hysteresis loop of AMa{H) which is defined 
by at Ma{H) at 3.3 K minus a diamagnetic background 
given by XdH with Xd — -7.317 x 10"'' emu/gOe. As 
H decreases, a trapped magnetic flux, corresponding to 
a paramagnetic moment Mr (« 2.7 x 10^^ emu/g) re- 
mains in the sample. With further cycling of H from -500 
to 500 Oe a characteristic hysteresis loop is observed. 

We also measured the hysteresis loop of the magneti- 
zation Mc{H) for H\\c at T = 1.9 K. The difference AM^ 
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FIG. 4: Hysteresis loop of DC magnetization minus a dia- 
magnetic background (= XdH) for Bi-MG. H±c. (a) T = 1.9 
K. Xd = -5.143 X lO"'^ emu/gOe. (b) T = 3.3 K. Xd = -7.317 
x 10~^ emu/gOe. 



is defined by Mc{H) at 1.9 K minus a diamagnetic back- 
ground given by XdH with Xd = -2.980 x 10~® emu/gOe. 
The hysteresis loop of AMc for H\\c is very similar to 
the corresponding data AMa for HJ-c. The value of Hd 
{— 75 Oe) for H\\c is a little lower than that for HLc. 
The peak value of AMc is almost the same as that of 
AMa- The hysteresis loop of AMc at 1.9 K exhibits a 
nearly reversible behavior. In bulk samples this could be 
explained in terms of a lack of structural defects to pro- 
vide pinning sites in the vortex state {Hd < H < Hc2)- 
In Bi-MG, however, the size of islands is much shorter 
than the coherence length ^. Thus the pinning effect is 
not relevant for these islands. In contrast, the hysteresis 
loop of AMa at 1.9 K exhibits some hysteretic behavior. 
Note that similar magnetization curves are observed in a 
superconductor TaC nanoparticles which are synthesized 
using a vapor-solid reaction path starting with carbon 
nanotube precursor^ 



B. Curie-like susceptibility due to nanographites 

Figure [5l shows typical data of xfc for H = 1 kOe for 
H±c and H\\c. The susceptibility is negative at high T 
and drastically increases with decreasing T. The least 
squares fit of the data {xfc vs T) for 1.9 < T < 30 K to 
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FIG. 5: T dependence of FC susceptibilities for Bi-MG. H 
— 1 kOe. HJ-C and H\\c. The inset shows the reciprocal 
susceptibility {x — Xo)~^ as a function of T. The solid lines 
are described by Eq.(|2J with parameters given in the text. 
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FIG. 6: T depedence of c-axis resistivity pc{T,H) for Bi-MG 
at various H. H\\c. I\\H. 



the Curie- Weiss law 

C 

X^Xo+ y _Q , (2) 

yields 9 = -0.86 ± 0.05 K, C = (1.999 ± 0.029) x 10"^ 
emu K/g, xo = (-1.307 ± 0.001) x 10"^ emu/g for H±c, 
and 6 = -1.13 ± 0.02 K, C = (1.460 ± 0.008) x lO''^ emu 
K/g, Xo = (-1.146 ± 0.001) x 10"^ emu/g for H\\c. The 



value of Q is very close to zero, showing a Curie- like law. 
In the inset of Fig. \^we show the reciprocal susceptibil- 
ity (x ^ Xo)^^ a fmiction of T. The negative value xo 
is from the orbital diamagnetic susceptibility. There is 
a crossover from the high-temperature diamagnetic sus- 
ceptibility to the low-temperature Curie-like susceptibil- 
ity around 50 K. We assume that the susceptibility of 
Bi-MG at 100 K corresponds to the diamagnetic suscep- 
tibility since the Pauli susceptibility is positive and is 
nearly equal to zero at 100 K. From Fig. El the diamag- 
netic susceptibility for H\\c and i/_Lc is estimated as Xc ~ 
-1.1 X 10~^ emu/g and Xa ~ - 1-3 x 10~^ emu/g, which 
are almost isotropic. These values are in contrast with 
the susceptibility of HOPG, which is very anisotropic: Xc 
= -25.86 X 10~^ emu/g and Xa = -1.06 x 10"'^ emu/g at 
H = 1 kOe. The absolute value of Xc in Bi-MG is much 
smaller than that in HOPG, while the values of Xa are 
on the same order for both systems. 

Here it is interesting to compare our data of suscepti- 
bility with that of nanographites prepared by heat treat- 
ing diamond nanoparticle^^ and a disorder network of 
nanographites in activated carbon fiber (ACF).^ The T 
dependence of the susceptibility for these compounds is 
similar to that of Bi-MG. In particular, the values of 0, 
C, and Xo for the ACF are on the same order as those 
for Bi-MG, where 9 « -2 K, C = 1.21 x lO"'^ emu K/g, 
and Xo = -1-36 x 10~^ emu/g for ACF prepared at the 
heat treatment temperature HTT = 1500 °C. — Here we 
assume that the spin S" is 1/2 and the Lande g-factor 
is 2 since the spins are associated with carbon materials. 
The number of spin density per 1 g of Bi-MG is estimated 
as (4.7 - 6.4) x 10^^/g, which is comparable with 3.9 x 
lO^Vg for the ACF with HTT = 1500 °Gi^ 

It has been theoretically predicted by Wakabayashi et 
al}^ that the Pauli susceptibility exhibits a Curie- like be- 
havior in the nanographites with zigzag edges, because 
of the DOS which is sharply peaked at the Fermi en- 
ergy. However, qualitatively we think that the Curie-like 
behavior at low T in Bi-MG is due to the conduction elec- 
trons localized around the zigzag edges of nanographites, 
which have local magnetic moments (spin S* = 1/2 and 
a Lande g-factor g = 2). The origin of spin polarization 
in nanographites with zigzag edges has been discussed 
by Fujita et al.^'^ using the Hamiltonian that consists 
of the on-site Coulomb repulsive interaction (U) when 
the site is occupied by two electrons, and the electron 
transfer integral between the nearest sites (t) . They have 
shown that a ferrimagnetic spin polarization emerges on 
the edge carbons even for weak U/t « 0.1. 

The orbital diamagnetic susceptibility is very sensitive 
to the size and edge shapes of nanographites. The orbital 
diamagnetic susceptibility is almost isotropic in Bi-MG, 
while it is very anisotropic in pristine graphite with infi- 
nite size. In Bi-MG the orbital cyclotron motion of elec- 
trons in the presence of H {H\\c) is greatly suppressed. 
This result is in good agreement with the prediction 
by Wakabayashi et alii^ In ribbon-shaped nanographites 
with zigzag edges, the magnitude of the diamagnetic sus- 
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FIG. 7: T dependence of in-plane resistivity pa{T,H) for Bi- 
MG at // = 0, 100, 200, and 400 Oe. H±c. I±H. The inset 
shows the T dependence of Pa{T,H = 0) (•) and dpa{T, H 
= 0)/dr (O) at low T. 

ceptibility decreases as the ribbon width decreases. The 
flow of the orbital diamagnetic ring current significantly 
depends on the lattice topology near the graphite edge. 



C. c-axis and in-plane resistivities 

The c-axis resistivity pc of Bi-MG was measured us- 
ing the four-probe method. After the sample was cooled 
from 298 K to 1.9 K at = 0, the c-axis resistivity 
Pc was measured with increasing T from 1.9 K to 50 K 
without and with H [H\\c (c: c axis)], where the current 
direction is parallel to the field direction. This resistivity 
is denoted as the longitudinal magnetoresistance. Figure 
El shows the T dependence of pc for Bi-MG a,t H — - 
10 kOe. The resistivity pc for Bi-MG shows a very broad 
peak around 14 K at 7? = 0, in contrast to pc for HOPG 
exhibiting a peak around 40 Kii^ This peak shifts to the 
low T side with increasing H. For H > 20 kOe, pc shows 
a semiconductor-like behavior: it decreases with increas- 
ing T. Note that pc shows a positive magnetoresistance: 
it increases with increasing H at any fixed T between 1.9 
and 50 K. The value of pc at T = 298 K and 7? = is 0.11 
57cm, which is on the same order as that of HOPG {pc ~ 
0.096 ncm)i^ and BiClg GIG {pc = 0.2 f^cm)," leading 
to a mean free path less than 1 A according to the Drude 
formula. This result suggests that there is no overlapping 
over nearest-neighbor layers along the c axis. The c-axis 
conduction can occur via the hopping of carriers between 
layers. 

The in-plane resistivity pa of Bi-MG was also mea- 
sured using the four-probe method. After the sample 



was cooled from 298 K to 1.9 K at 77 = 0, the in-plane 
resistivity pa was measured with increasing T from 1.9 K 
to 50 K without and with H {H±c), where the field di- 
rection is perpendicular to the current direction in the c 
plane. This resistivity is denoted as the transverse mag- 
netoresistance. The resistivity ratio A Pc/ Pa) is es- 
timated as 30 at 298 K using the measured pa- How- 
ever, the actual value of A is considered to be much 
larger than 30 because of possible contribution of pc to 
the measured pa- Figure [7| shows the T dependence of 
Pa of Bi-MG at 7? = - 400 Oe. The zero-field [H 
= 0) in-plane resistivity increases with increasing T at 
low T, showing a metallic behavior. It has a maximum 
around 15 K, and it decreases with further increasing 
T, showing a semiconductor-like behavior. Note that 
the value of pa for Bi-MG 4.2 milcm) is much larger 
than that of BiGlg GIG (= 27 p^cm at 298 K). We find 
that the T dependence of pa for Bi-MG is very similar 
to that of grafoil: they even have similar magnitudes. 
The grafoil is a pyrolytic graphite with a polycrystalline 
structures formed of many domains. According to Koike 
et al.)2i Pa of grafoil (grade GTA, Union carbide) shows 
a semiconductor-like behavior, while pa of HOPG (Union 
Garbide) shows a metallic behavior. This result sug- 
gests that the semiconductor-like behavior in Bi-MG may 
arise from nanographites where the degree of disorder is 
greatly enhanced. Such a semiconductor-like behavior is 
observed at least below 25 kOe. 

In the inset of Fig. [7| we show the detail of pa at H 
= near Tc = 2.48 K. No drastic decrease in pa below 
Tc is observed with decreasing T, while the T-derivative 
dpa/dT gradually decreases with increasing T around Tc. 
The causes for the finite value of pa below Tc in spite of 
the superconducting phase are as follows, (i) The sample 
used in the present work is an exfoliated Bi-MG based 
on HOPG. The sample surface is not flat partly because 
of cracks generated in the basal plane. Since the current 
path is not always located on the same layer, the con- 
tribution of large pc to the observed pa is not negligibly 
small. Another possibility is the local superconductivity 
in isolated islands. In such a system there is competi- 
tion between the Josephson coupling (Ej) and the charg- 
ing energy (Eu) between superconducting islands «2i. For 
Ej ^ Ejj, the Gooper pairs are delocalized leading to 
a superconducting state with vanishing resistivity. For 
Ej <^ Eu, the pairs will be localized and the transport 
is possible only by thermal activation leading to insulat- 
ing behavior at T = K. The resistivity varies with T 
as p ~ eyi'p{Ta/T), where Ta is related to the activation 
energy. In this model, the resistivity should increase with 
decreasing T, which contradicts with our result. Liu et 
al.^^ have reported an unusual T dependence of resistiv- 
ity at low H in ultrathin superconducting films of Pb, Al, 
and Bi. The resistivity varies with T as p ~ exp(r/To) 
at low r, where Tq is a characteristic temperature. In Bi- 
MG, dp/dT decreases with increasing T around T^- In 
superconducting thin films, dp/dT exp(T/ro)/ro] in- 
creases with increasing T . Therefore, the T dependence 
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FIG. 8: H dependence of pa{T,H) for Bi-MG at fixed T (2 
< r < 50 K). H±c. 



of Pq in Bi-MG is not the case of p « exp(r/ro). 



D. 2D Weak localization effect 

As shown in Fig. [T] the in-plane resistivity pa sUghtly 
decreases with increasing H at the same T (at least for 
1.9 < T < 50 K), indicating a negative magentoresistance 
(NMR). Figures IHIa) and (b) show the H dependence of 
Pa for Bi-MG with 2 < T < 50 K, where H±c and the 
field direction is perpendicular to the current direction. 
For each T, pa decreases with increasing H at low H 
(H > 100 Oe), exhibits a local minimum around H — 
2.5 kOe, and increases with further increasing H. The 
sign of the difference Apa [= Pa{T, H) - pa{T, H = 0)] is 
negative for < H < 3.5 kOe due to the possible 2D weak 
localization effect (WLE)i^ It changes from negative to 
positive at i? = 3.5 kOe. The resistivity pa drastically in- 
creases with further increasing iJ, as is observed in com- 
pensated metals such as bulk Bip^ The H dependence of 
Pa for 3 kOe < H < A5 kOe can be well described by pa 
= po + pi-ff + P2H^, where po, pi, and p2 are constants: 
po = 4.290 mf^cm, pi = (9.85 ± 0.11) x lO'^ /if]cm/Oe, 
and p2 = (4.56 ± 0.23) x 10"^ pQcm/Oe^ at 1.9 K. The 
linear term {piH) is dominant compared to the squared- 
power term [piH"^) for H <^ Hi, where Hi = pi/p2 — 
216 kOe at 1.9 K. For comparison, we also measured the 
H dependence of pa in Bi-MG for 1.9 < T < 30 K, where 
H\\c. At each T, pa increases with increasing H. Thus 



the sign of Apa is always positive for < H < 47.5 kOe 
and 1.9 < r < 50 K, indicating an anti-localization effect 
which is similar to that observed in Bi thin films^^ with 
strong spin-orbit interactions. 

The T dependence of pa at H = 45 kOe for H±c is 
described by 



Pa ao - ai ln(T), 



(3) 



for 1.9 < T < 4.3 K, where oq = (4.8261 ± 0.0002) mf^cm 
and ai = (3.79 ± 0.19) x 10^^ mf2cm. As shown in 
previous papers^iSiS the theory of the 2D WLE predicts 
that the following relation is valid for the ratio a2/ai, 



ai 



(4) 



with A = ap+-f, where e^/{2T:'^fi) = 1.23314 x 10"^ 
and is the in-plane conductivity defined by Id p\ {p'i, 
is the in-plane resistivity and Ic is the c axis repeat dis- 
tance). In the parameter A, a is nearly equal to unity, 
p is the exponent of the inelastic life time of the con- 
duction electrons (r^ w T~p), and 7 is a numerical factor 
giving a measure of the screening by other carriers. For 
convenience, here we use the value of ct^d (= 1-89 x 10^^ 
il^^) for kish graphite {Ic — 3.35 A and p" = 1.77 pflcm 
at T = 4.2 K) obtained by Koike et al.yS.'^ instead of the 
corresponding data for Bi-MG. A kish graphite is a high- 
quality single crystal of graphite which is deposited on 
walls of blast furnaces for steel production. Using ai/ao 
— 7.853 X 10~^ the parameter A is estimated a.s A — 1.20, 
which is comparable with 1.14 for stage-4 M0CI5 GICi^ 
These results suggest that the logarithmic behavior of pa 
can be explained in terms of the 2D WLE. 



E. Field induced metal-semiconductor transition 

In Fig. inia) we show the T dependence of the differ- 
ence (pa) between pa at T and that at 20 K for various 
H (Hlc): Ha = PaiT,H) - PaiT ^ 20K,H). The T 
dependence of p,a below 10 K is dependent on H. The 
difference pa shows a metallic behavior (dpa/dT > 0) 
below 15 - 20 kOe, while it shows a semiconductor-like 
behavior {dpa/dT < 0) above H ^ 25 kOe. In Fig. E^b) 
we show the T dependence of pa for H\\c. The T depen- 
dence of Pa below 7 K is dependent on H. The difference 
Pa shows a metallic behavior {dpa/dT > 0) below 10 
kOe, while it shows a semiconductor-like behavior above 
20 kOe. Note that pa ai H = 47.5 kOe drastically in- 
creases with decreasing T, leading to the insulating state. 
The magnitude of pa for iJ ||c is on the same order as that 
for HLc. The T dependence of pa at H = 47.5 kOe for 
H\\c is well described by Eq.® for 1.9 < T < 4.3 K, 
where oq = (5.8385 ± 0.0002) mflcm and oi = (5.33 ± 
0.23) X 10~^ mficm. The corresponding fitting curve is 
denoted by the solid line in the inset of Fig. EJb). The 
ratio fli/ao = 9.13 x 10""* is almost the same as that (= 
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7 85 X 10-'') for paSX H = 45 kOe for ff_Lc (see Sec. 
1111 C|) . indicating the occurrence of 2D WLE. 

Bi-MG undergoes a transition from the metaUic phase 
to the semiconductor-hke phase at a critical field He- 
Similar crossover behavior is observed in Si MOSFET)^ 
and amorphous metal films of InO^;;^ MoGe^ and Bii^ 
The value of (« 25 kOe) for Bi-MG is almost the 
same as that for Si MOSFET^ and amorphous MoGe 
filmii Note that the Zeemann energy gSfisH at 25 kOe 
corresponds to a thermal energy ksTn with Th — 1-7 K, 
where g = 2 and 5 = 1/2. The temperature Th is slightly 
lower than Tc (= 2.48 K) for the superconductivity. The 
suppression of the metallic phase by H is independent of 
the directions of H {H\\c and H±c) for Bi-MG. These 
results suggest that the spin related effect is significant 
compared to the orbital effect. This is consistent with 
the result derived from the susceptibility measurement 
that the Curie-like susceptibility is dominant at low T. 



Bi-MG 




FIG. 9: T dependence of the difference fia between Pa{T, H) 
and pa{T = 20K,H) [i^a = Pa{T,H) - pa{T = 20K,H)] for 
Bi-MG with various H {I < H < 45 kOe). ILH. (a) HLc. 
(b) H\\c. The inset of (b) shows the T dependence of pa{T, H) 
at H — 47.5 kOe. The solid line denotes the fitting curve of 
the data to Eq.J^J. 



IV. CONCLUSION 

Bi-MG shows superconductivity at Tc = 2.48 K, where 
the coherence length is much larger than the size of 
nanoparticles. The spin related effect characterized by 
a Curie-like susceptibility is enhanced, while the orbital 
effect is greatly suppressed. The transition from metal- 
lic phase to semiconductor-like phase is induced by the 
application of H above 25 kOe. A negative magnetore- 
sistance in pa and a logarithmic divergence of pa with 
decreasing T are indicative of the 2D WLE. 



Acknowledgments 

We would like to thank Kikuo Harigaya for valu- 
able comments on the magnetism of nanographites with 
zigzag edges. The work at SUNY-Binghamton was 
supported by SUNY-Research Foundation (Grant No. 
240-9522A). The work at Osaka University was sup- 
ported by the Ministry of Cultural Affairs, Education 
and Sport, Japan under the grant for young scientists 
(No. 70314375) and by Kansai Invention Center, Kyoto, 
Japan. We were grateful to Advanced Ceramics Corpo- 
ration, Ohio for providing us with HOPG (grade ZYA). 



* |suzuki@binghamton. edu| 

^ E. Abrahams, P.W. Anderson, D.C. Licciardello, and T.V. 

Ramakrishnan, Phys. Rev. Lett. 42, 673 (1979). 
^ S.V. Kravchenko, W.E. Mason, G.E. Bowker, J.E. 

Furneaux, V.M. Pudalov and M. D'lorio, Phys. Rev. B 

51, 7038 (1995). 
^ D. Simonian, S.V. Kravchenko, M.P. Sarachik, and V.M. 

Pudalov, Phys. Rev. Lett. 79, 2304 (1997). 

* S.V. Kravchenko and T.M. Klapwijk, Phys. Rev. Lett. 84, 
2909 (2000). 



^ A.F. Hebard and M.A. Paalanen, Phys. Rev. Lett. 65, 927 
(1990). 

^ Ah Yazdani and A. Kapitulnik, Phys. Rev. Lett. 74, 3037 
(1995). 

N. Mason and A. Kapitulnik, Phys. Rev. Lett. 82, 5341 
(1999). 

* N. Markovic, C. Christiansen, and A.M. Goldman, Phys. 

Rev. Lett. 81, 5217 (1998). 
^ J. Walter, H. Shioyama, and Y. Sawada, Carbon 36, 1811 

(1998). 



9 



1° J. Walter and H. Shioyama, Carbon 37, 1151 (1999). 
" E. McRao and J. Walter, Carbon 39, 717 (2001). 

M. Suzuki, C.Lee, I.S. Suzuki, K. Matsubara, and K. Sug- 

ihara, Phys. Rev. B 54, 17128 (1996). 

M. Pujita, K. Wakabayashi, K. Nakada, and K. Kusakabe, 

J. Phys. Soc. Jpn. 65, 1920 (1996). 

K. Wakabayashi, M. Sigrist, and M. Fujita, J. Phys. Soc. 
Jpn. 67, 2089 (1998). 

K. Wakabayashi, M. Fujita, H. Ajiki, and M. Sigrist, Phys. 
Rev. B 59, 8271 (1999). 

J. Walter and H. Shioyama, Phys. Lett. A 254, 65 (1999). 

J. Walter, Adv. Mater. 12, 31 (2000). 
1* J. Waher, Phil. Mag. Lett. 80, 257 (2000). 

D.B. Haviland, Y. Liu, and A.M. Goldman, Phys. Rev. 

Lett. 18, 2180 (1989). 
^° B. Weitzel and H. Micklitz, Phys. Rev. Lett. 66, 385 

(1991). 

G. Deutscher, O. Eutiu-Wohlman, and Y. Shapira, Phys. 
Rev. B 22, 4264 (1980). 

D. Saint-James and P.G. de Gennes, Phys. Lett. 7, 306 

(1963). 

S. Tanuma, in Graphite Intercalation Compounds 11, edited 
by H. Zabel and S.A. Sohn (Springer- Verlag, Berhn, 1992) 
p.l63. 



A. Fukunaga, S. Chu, and M.E. McHenry, J. Mater. Res. 

13, 2465 (1998). 

O.E. Andersson, B.L.V. Prasad, H. Sato, T. Enoki, Y. 
Hishiyama, Y. Kaburagi, M. Yoshikawa, and S. Bandow, 
Phys. Rev. B 58, 16387 (1998). 
^® Y. Shibayama, H. Sato, T. Enoki, and M. Endo, Phys. 
Rev. Lett. 84, 1744 (2000). 

Y. Koike, S. Morita, T. Nakamoto, and T. Fukase, J. Phys. 
Soc. Jpn. 54, 713 (1985). 

W. Zwerger, in Quantum Coherence in Mesoscopic Sys- 
tems, edited by B. Kramer (Plenum Press, New York, 
1991) p.491. 

Y. Liu, D.B. Haviland, L.L Glazman, and A.M. Goldman, 
Phys. Rev. Lett. 68. 2224 (1992). 

S. Hikami, A.L Larkin, and Y. Nagaoka, Prog. Theor. 
Phys. 63, 707 (1980). 

M. Suzuki and S. Tanuma, J. Phys. Soc. Jpn. 45, 1645 
(1978). 

F. Komori, S. Kobayashi, and W. Sasaki, J. Phys. Soc. 
Jpn. 52, 368 (1983). 

M. Suzuki, I.S. Suzuki, K. Matsubara, and K. Sugihara, 
Phys. Rev. B 61, 5013 (2000). 



